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Abstract. In this article, we analyze the microlocal properties of the linearized forward scattering operator 
F and the normal operator F*F (where F* is the adjoint of F) which arises in Synthetic Aperture Radar 
' imaging for the common midpoint acquisition geometry. When F* is applied to the scattered data, artifacts 

^Nj , appear. We show that F*F can be decomposed as a sum of four operators, each belonging to a class of 

distributions associated to two cleanly intersecting Lagrangians, /^''(Ao, Ai), thereby explaining the latter 
O ' artifacts. 

o 



< 



1. Introduction 



In this article, we analyze the microlocal properties of a transform that appears in Synthetic Aperture 
Radar (SAR) imaging. In SAR imaging, a region on the surface of the earth is illuminated by an 
electromagnetic transmitter and an image of the region is reconstructed based on the measurement of 
scattered waves at a receiver. For an in-depth treatment of SAR imaging, we refer the reader to [ll|2]. The 
CN ■ transform we study appears as a result of a common midpoint acquisition geometry: the transmitter and 
^ ! receiver move at equal speeds away from a common midpoint along a straight line. This geometry is of 
interest in bistatic imaging and in certain multiple scattering scenarios |14j . We first consider the linearized 
scattering operator F and show that it is an FIG. Since the conventional method of reconstructing the 
image of an object involves "backprojecting" the scattered data, we next study the composition of F with 
its adjoint F*. One of the main goals of this article is to understand the distribution class of F*F. 

In general the composition of two FIOs is not an FIG. Gne needs additional geometric conditions 
such as the transverse intersection condition |10j or the clean intersection condition [1] to make the 
composition operator again an FIG. When these assumptions fail to be satisfied, it is very useful to 
study the canonical relation associated to an FIG by considering the left and the right projections. More 
^ ■ precisely, let X and Y be manifolds and let Y; C) be the class of Fourier integral operators (FIGs) 

F : £'{X) — )• T)'{Y) of order m associated to the canonical relation C C {T*Y x T*X) \ {0} and denote by 
: C — 7- T*y, ttr : C — )• T*X, the left and right projections respectively. Where and how these projections 
drop rank determine the nature of the normal operator F*F. 

The singularities which appear in previous work related to SAR \13\ [5l [H [11] are folds and blowdowns, 
that is, ttl and vr/j have both fold singularities or ttl has a fold singularity and vr/j has a blowdown 
singularity. These singularities will be defined in Section [3l Then it is known that the corresponding 
normal operator belongs to a class of distributions /^™'''(A,C') introduced by [9] (and defined in Section 
[3]). This means that the adjoint operator F* introduces an additional singularity given by C apart from 
the initial one given by A. In this article, the linearized scattering operator F exhibits a new feature: both 
projections drop rank by one on a disjoint union of two smooth hypersurfaces Si U S2. Gn each of them, 
TTi is a projection with fold singularities and vr/j is a projection with blowdown singularities. We then show 
that F*F belongs to the class /2'"'0(A, Ci) + P'"'^{A,C2) + /2™'0(C7i, C3) + ^"^'^{02,03) (where these 
classes are given in Definition 13. 6p . This means that the adjoint operator F* adds three more singularities 
given by Ci, C2, C3 in addition to the true reconstructed singularity given by A. We clarify this in detail in 
Section [5l The main tool for proving our result is the iterated regularity property; a characteristic property 
of JP'^ classes [TJ Proposition 1.35]. 
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2. Statement of the main results 



2.1. The linearized scattering model. For simplicity, we assume tliat botii tlie transmitter and receiver 
are at the same height h > above tlie ground at all times and move in opposite directions at equal 
speeds along the line parallel to the x axis from the common midpoint (0, 0,/i). Such a model arises 
when considering signals which have scattered from a wall within the vicinity of a scatterer and can be 
understood in the context of the method of images; see [13] for more details. 

Let 7t(s) = {s,0,h) and 7_r(s) = (— s,0, /i) for s £ (0,oo) be the trajectories of the transmitter and 
receiver respectively. 

The linearized model for the scattered signal we will use in this article is 



d{s,t) := FV{s,t) = I e 'o^^''''^'> a{s,x,uj)V{x)dxduj for {s,t) G (0,oo) x (0,oo), (1) 



where V{x) = V{xi,X2) is the function modeling the object on the ground, R{s,x) is the bistatic distance: 

R{s,x) = |7t(s) - x\ + \x- 7ij(s)|, 
Co is the speed of electromagnetic wave in free-space and the amplitude term a is given by 

a{s,x,uj) = — — — — -, (2) 

167r^|7T(s) - x\\^r{s) - x\ 

where p is the Fourier transform of the transmitted waveform. 

2.2. Preliminary modifications on the scattered data. For simplicity, from now on we will assume 
that Co = 1. To make the composition of F with its adjoint F* to be well-defined, we multiply d{s,t) 
by an infinitely differentiable function f{s,t) identically 1 in a compact subset of (0,oo) x (0, oo) and 
supported in a slightly bigger compact subset of (0, oo) x (0, oo). We rename f ■ d as d again. 

As we will see below, our method cannot image a neighborhood of the common midpoint. That is, if 
the transmitter and receiver are at (s, 0, h) and (— s, 0, h) respectively, we cannot image a neighborhood of 
the origin on the plane. Therefore we modify d further by considering a smooth function g{s,t) such that 



g{s, t) = for (s, t) : \t - 2 vV + Zi^l < 20e^/h, (3) 

where e > is given. Again we let g ■ d to be d and g ■ a to he a. The choice of constant on the right hand 
side of ^ will be justified in Appendix iBl 

d{s, t) = FVis, t) = J ) ^(,^ ^^vix) dxd0 (4) 

and the phase function is 

(p{s, t, x,uj) = uj ^ — {xi — sY + X2 + h? — ^ (xi -I- sY + X2 + h"^ ^ . (5) 

From now on, we will denote the {s,t) space as Y and x = {xi,X2) space as X. 

We assume that the amplitude function a G S"''^^, that is, it satisfies the following estimate: For every 
compact set K C Y x X , non-negative integer a, and 2-indexes /3 = (/3i,/32) and 7, there is a constant c 
such that 

\dy^'dtd2a{s,t,x,u:)\ < c(l + |^|)-+(i/2)-. (6) 

This assumption is satisfied if the transmitted waveform from the antenna is approximately a Dirac delta 
distribution. 

With these modifications, we show that F is a Fourier integral operator of order m and study the 
properties of the natural projection maps from the canonical relation of F. Our first main result is the 
following: 

Theorem 2.1. Let F be as in @. Then 
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(a) F is an FIO of order m. 

(h) The canonical relation C associated to F is given by 

Q ^ i ( g ^ ^1 - ■s xi + s \ 

\ \ ' ' ^y'ixi- sy + x^ + h^ Y^(xi + s)2 + x^ + /i2/' 

/ Xl - S Xl + s \ 

xt,X2,—uj\ — 1 (7) 

V^(xi -s)2 + x2 + /l2 ^(xi + s)2 + x2 + /l2;' ^ 

— a;( — ^ = H ^ ^ ) I where 

^y/{xi- s)'^ + xl + h'^ V(xi+s)2 + x| + /l2/ / 



s > 0, t = y' (xi - s)2 + + /i2 + y (xi + s)2 + x| + /i2, X / 0, and / 0| , 

and C /las global parameterization (0, oo) x (M2 \ O) x (M \ 0) 9 (s, xi, X2, w) i— )• C. 

(c) Let ttl '■ C ^ j'*Y and iipi, : C — >■ T*X be the left and right projections respectively. Then 
ttl and ttr drop rank simply by one on a set S = Si U S2 where in the coordinates {s,x,u}), 
Si = {(s, Xl, 0, ti;)|s > 0, |xi| > e', w / 0} and S2 = {(s, 0, X2, > 0, |x2| > e' ,u 7^ 0} for < e' 
small enough. 

(d) ttl has a fold singularity along S. 

(e) TTji has a blowdown singularity along S. 

Remark 2.2. Note that due to the function g{s,t) of ^ m the amphtude, it is enough to consider only 
points in C that are strictly away from {{s,0,uj) : s > 0,uj ^ 0}. This is reflected in the definitions of Si 
and S2, where |xi| and |x2|, respectively, are strictly positive. 

Remark 2.3. Note that C is even with respect to both xi and X2. In other words C is a 4-1 relation. This 
observation suggests that ttl (respectively ttjj) has two fold (respectively blowdown) sets. See Proposition 

We then analyze the normal operator F*F. Our next main result is the following: 

Theorem 2.4. Let F be as in Q of order m. Then F*F can be decomposed into a sum belonging to 
/2'".0(A,Ci) + /2™'°(A,C2) + /2"''°(Ci,C3) + l2™'°(C2,C3) where these classes are given in Definition\EE 

In Remark 15.71 we will explain why the added singularities given by Ci , C2 , C3 have the same strength 
as the object singularities given by A. 

3. Preliminaries 



3.1. Singularities and classes. In this section we will define fold and blowdown singularities and 
describe the class of distributions required for the analysis of the composition operator F*F. 

Definition 3.1 ([H p. 109-111]). Let M and be manifolds of dimension n and let / : M be C°°. 

Let be a non- vanishing volume form on N and define T, = {a € M : f*Q(a) = 0}, that is, S is the 
set of critical points of /. Note that, equivalently, S is defined by the vanishing of the determinant of the 
Jacobian of /. 

(a) If for all cr G S, we have (i) the corank of / at u is 1, (ii) ker(d/o-) n T^-S = {0}, (iii) vanishes 
exactly to first order on S, then we say that / is a fold. 

(b) If for all fj G S, we have (i) the rank of / is constant; let us call this constant k, (ii) ker(d/o-) C Tq-S, 
(iii) f*Q vanishes exactly to order n — A; on S, then we say that / is a blowdown. 

We now define /^'^ classes. They were first introduced by Melrose and Uhlmann jl2j . Guillemin and 
Uhlmann [9] and Greenleaf and Uhlmann ^ and they were used in the context of radar imaging in |13t El |6] . 
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Definition 3.2. Two submanifolds M and N intersect cleanly if M n iV is a smooth submanifold and 
T{M r\N) = TM n TN. 

Let us consider tlie following example 

Example 3.3. Let Aq = At-R" = {(x, ^; x, e M", ^ e \ 0} he the diagonal in T*W x T*W and let 
Ai = {{x' ,Xn,^' ,^]x' ,ymi' ,Q)\x' G M""^, ^' G \ 0}. Then, Aq intersects Ai cleanly in codimension 

1. 

Now we define the class of product-type symbols S'^'\m,n,k). 

Definition 3.4. SP'\m, n, k) is the set of aU functions a{z, ^, a) G C°°(]R™ x M" x M.'') such that for every 
K C and every a G /3 G 7 E there is ck a s y such that 

|a,"af5Ja(z,e,a)| < CK,a,is,y{l + + k|)'-l^l,V(z,C,T) G X X 

Since any two sets of cleanly intersecting Lagrangians are equivalent [9], we first define I^'' classes for 
the case in Example 13.31 

Definition 3.5. [9] Let /^''(Aq, Ai) be the set of all distributions u such that u = U1+U2 with ui G 
and 

with a G 5^'''' where p' = p - f + | and I' = I - \. 

Let (Ao,Ai) be a pair of cleanly intersection Lagrangians in codimension 1 and let x t>e a canonical 
transformation which maps (Aq, Ai) into (Aq, Ai) and maps AoH Ai to AoHAi, where A^ are from Example 
Next we define the /^''(Ao, Ai). 

Definition 3.6 ([2]). Let /^''(Ao, Ai) be the set of all distributions u such that u = ui + U2 + '^Vi where 
ui G /^+'(Ao \ Ai), U2 G /^(Ai \ Aq), the sum ^i^j is locally finite and Vi = Awi where A is a zero order 
FIO associated to x^^j the canonical transformation from above, and Wi G /^''(Aq, Ai). 
If u is the Schwartz kernel of the linear operator F, then we say F G /^''(Aq, Ai). 

This class of distributions is invariant under FIOs associated to canonical transformations which map 
the pair (Aq, Ai) to itself. If F G /P''(Ao, Ai) then F G /p+'(Ao \ Ai) and F G /^(Ai \ Aq) [9]. Here by 
F G /^"'"'(Aq \ Ai), we mean that the Schwartz kernel of F belongs to I^"*"^ ( Aq \ Ai ) microlocally away from 
Ai. 

One way to show that a distribution belongs to class is by using the iterated regularity property: 

Proposition 3.7. [3 Proposition 1.35] If u G V'{X x Y) then u G /^''(Ao, Ai) if there is an sq £ R such 
that for all first order pseudodifferential operators Pi with principal symbols vanishing on Aq U Ai, we have 

P,P2...PrU£Ht:^^. 



4. Analysis of the Operator F 
In this Section, we prove Theorem l2.lt as a result of Lemma l4.1l and Proposition 14.31 
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Lemma 4.1. F is an FIO of order m with the canonical relation C given by 

— u}{—= -I = ) I wit/i 

-s)2 + x| + /l2 y/(xi +s)2 +2;2 + /l2y / 



s > 0,f = - s)2 + xi + /i2 + ^{xi + s)2 + xi + /i2,x G \{0},tj / 0|. 

We note that (0, oo) x (M2 \^ q) x (M \ 0) 9 (s, xi, X2, w) ^ C is a global parametrization of C . 

We will use the coordinates {s^x^uj) in this lemma from now on to describe C and subsets of C. 

Proof. The phase function ip is non-degenerate with dx^, dg^t^ nowhere whenever duj^p = 0. We should 
mention that Vdco^p ^ 0. (Note that in order for dx^ to be nowhere 0, we require exclusion of the common 
midpoint from our analysis). This observation is needed to show F is a FIO rather than just a Fourier 
integral distribution. Recalling that a satisfies amplitude estimates ([6|), we conclude that F is an FIO |15j . 
Also since a is of order m + ^, the order of the FIO is m [3l Definition 3.2.2]. By definition [101 Equation 
(3.1.2)] 

C = {{{s,t,ds(p,dt(p);{x,-dx(p)) : dujip = 0}. 
A calculation using this definition establishes dH). Furthermore, it is easy to see that (s, xi, X2, w) is a 
global parametrization of A. □ 

Remark 4.2. In the SAR application, a has order 2 which makes operator F of order |. But from now 
on will consider that F has order m. 

Proposition 4.3. Denoting the restriction of the left and right projections to C by ttl andnn respectively, 
we have 

(a) TTL and ttr drop rank by one on a set S = Si U S2. Here we use the global coordinates from Lemma 

(b) TTL has a fold singularity along S. 

(c) TTji has a blowdown singularity along S. 

Proof Let A = y^(xi - s)2 + xl + h^ and B = a/(xi + s)2 + Xg + h"^. We have 

TrL{xi,X2, s, w) = {s,A + B, -( ^^^ ^ - — !-^^)w, -u) 

and 

/ 



and 



XI— S I Xl+S X2 I X2 

4 J B A ^ B 
,( xl+h? x^+h?-. ,(xi-s)x2 {x-i_+s)x2 ■ 
53 33 g3 , 





4xiX2sa;^^ , {xl — s"^ + x^ + h"^ . 



1 




* 





* 


* 





-1/ 



detdvr, = -i^(l+ - ^ ^i? ) (9) 

We have that s > and the number in the parenthesis is a positive number by Lemma 14.41 below. 

Therefore, vrx, drops rank by one on S = Si U S2. To show (i(det(d7rj^)) is nowhere zero on S, one uses 
the product rule in ([9]) and the fact that the differential of '^^^^^i^ is never zero on S and the inequality in 
Lemma [ 
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On Si the kernel of (Ittl is which is transversal to Si and on S2 the kernel of diri is which is 
transversal to S2. This means that ttl has a fold singularity along S. 
Similarly, 

Xl - S Xl+ S ,X2 X2. . 

vrii(xi,X2,s,a;) = (xi, X2, -(— \ ^)^'~(^ + 'b'^'' 



Then 



/l 

1 

* * -(^-^) -(^ + ^: 

V* * -..((-9)- - -(f + f) y 



has the same determinant so ttr drops rank by one on S and the kernel of dvr/j is a linear combination of 
-g^ and ^ which are tangent to both Si and S2. This means that ttr has a blowdown singularity along S. 

□ 

Lemma 4.4. For all s / 0, 

^ ^ xf - + xj + h'^ ^ ^ 
\x - 1t{s)\\x - -iR{s)\ 

Proof. Equivalently, we show that {\x — ^t{s)\\x — 7_r(s)|)^ > {x\ + x\ + h? — s^)^. Expanding out both 
sides and simplifying, we obtain 4s^(x2 + h?) > which holds for s 7^ 0, since h > 0. Therefore the lemma 
is proved. □ 



5. Analysis of the normal operator F*F 



We have 



X a{s, t, x, uj)a{s, t, y, uj)V {y)dsdtdujdujdy . 
After an application of the method of stationary phase in t and oj, the Schwartz kernel of this operator is 

K{x,y)= [ e''^(l"(^)l+l^~^«(^)l-(l^"^^(")l+l^-^«WI»5(x,y,s,w)dsda;. (10) 



Note that a G ^s^^+i since we assume a G 5""+i/2_ 

Let the phase function of the kernel K be denoted by 

^ = u{\y- 7t(s)| + \y- 7fl(s)| - i\x - 7t(s)| + \x - 7ij(s)|)) • (H) 

Proposition 5.1. The wavefront set of the kernel K of F*F satisfies, 

WF{K)' C A U C7i U C72 U C3, 

where A is the diagonal in T*X x T*X and the Lagrangians Ci for i = 1, 2, 3 are the graphs of the following 
functions Xi for i = 1, 2, 3 on T*X: 

Xi{x,0 = {xi,-X2,^i,-C2),X2{x,0 = (-2;i,2;2,-6,6) andx3 = Xi ° X2- 
Furthermore we have: 

(a) A and Ci, A and C2, Ci and C3, C2 and C3 intersect cleanly in codimension 2. 

(b) AnC3 = CinC2 = ^ 
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Proof. In order to find the wavefront set of the kernel K, we consider the canonical relation C* o C of F*F: 
C* o C = {(x, ^; y, r])\{x, ^; s, t, a, r) G C*; (s, t, cr, r; y, r/) G C}. We have that (s, t, a, r; y,r]) G C implies 



cr = r 



and (x, ^; s, cr, r) G C implies 



(yi - s) 


|2 + y2 + /j2 _^ ^ 




yi - s 




- sy + yl + h'^ 


-( 


yi- s 


V(yi 


- s)^ + y2 + h? 


-( 


y2 


V(yi 





yi + s 



Viyi + s)2 + + /i 
+ + ' ) (12) 

V(yi + «)' + y2' + ' 



t = (xi - s)2 + x| + /i2 + y + + x^ + /jS 
Xi — S Xi + s 



(T = T 



V^(xi - s)2 + + /l2 Y^(xi + S)2 + X^ + /l2 • 



( 2^1 — S ^ Xi + S ^ 

■\/(xi - s)2 + x| + /l2 Y^(xi + s)2 + x| + /l2 

6 = ^( , + — ) 

V^(xi - s)2 + x^ + /l2 y^(xi + S)2 + X^ + /l2 



From the first two relations in (1121) and (I13D. we have 



/(yi - + + + ^(yi + s)2 + y2 + = ^(^^ _ 5)2 + ^2 + /j2 + ^ (^^ + 5)2 + ^2 + /i2 (14) 

and 

yi- s yi + s xi- s xi + g ^^^^ 

Viyi - + yi + h'^ ^J{yi + sf + yl + h? ^{xi- sf + xl + h? ^'{xi + sf + xl + h'^' 
We will use the prolate spheroidal coordinates to solve for x and y. We let 

xi = s cosh p cos 4> yi = s cosh p' cos cp' 

X2 = s sinh p sin <j) cos ^ y2 = s sinh sin cos 9' (16) 
X3 = /i + s sinh p sin (j) sin ys = h + s sinh p' sin 0' sin 9' 

with p > 0, < (/> < TT and < 6* < 27r. 

In this case X3 = and we use it to solve for h. Hence 

(xi — s)^ + X2 + /i^ = s^(coshp — cos (/))^ 

and 

(xi + s)^ + X2 + /i^ = s^(coshp + cos 0)^ 
Noting that s > and cosh p it cos > 0, the first relation given by (|14p in these coordinates becomes 

s(cosh p — cos (p) + s(cosh p + cos cp) = s(cosh p' — cos (p') + s(cosh p' + cos </>') 
from which we get 

cosh p = cosh p' ^ p = p' 
The second relation given by psp becomes 

cosh p cos — 1 cosh p cos + 1 cosh p cos (p' — 1 cosh p cos + 1 
cosh p — cos (p cosh p + cos (p cosh p — cos (p' cosh p + cos cp' 
After simplification we get 

sin^ (p s\v? 6' 



cosh p — cos^ cosh p — cos^ 0' 
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which imphes 

(cosh^ p — l)(sin^ (f) — s\v? <^') = 

Thus sin (/> = it sin (p' ^ (j) = ±^', tt ± 

We remark that cos^ = ifcy^l — g2^i^^p^^^i ^ — =tcos0' and note that X3 = imphes that sin((^) 7^ 0, 

so that division by sin(i;^) is ahowed here. We also remark that it is enough to consider cos^ = cos0' as no 
additional relations are introduced by considering cos 9 = — cos 6'. 
Now we go back to x and y coordinates. 

If (f)' = (j) then xi = yi, X2 = 1/2, S,i = Vi foi' ^ = 1)2. For these points, the composition, 

C'oCcA = {{x,^-x,0}- 

If (j)' = —(j) then xi = Hi, —X2 = y2,(,i = —£.2 = ??2- For these points, the composition, C* o C is a 
subset of Ci = {(xi,X2,^i,C2;2;i, -X2,Ci, -6)} which is the graph of Xii^^O = i^i, -3^2,6) -'^2)- This in 
the base space represents the reflection about the xi axis. 

li (j)' = TT — (p then —xi = yi, X2 = y2, — = f?i,C2 = i]2- For these points, the composition C* o C is 
a subset of C2 = {(a^i, 2:2, Ci, 6; -xi,X2, -(.1,^2)} which is the graph of X2{x,£,) = {-xi,X2, -^1)6)- This 
in the base space represents the reflection about the X2 axis. 

If (p' = n + then —xi = yi, —X2 = y2, —^1 = "Hii = ??2- For these points, the composition, C*oC is a 
subset of C3 = {(xi,X2,6,6; -xi, -X2, -6, -6)} which is the graph of X3ix,0 = {-xi, -X2, -^i, -6)- 
This in the base space represents the reflection about the origin. 

Notice that xi o xi = Id, ^20X2 = Id, Xi°X2 = X3- 

So far we have obtained that C* o C C A U Ci U C2 U C3. 

Next we consider the intersections of any two of these Lagrangians. We have 

A intersects Ci cleanly in codimension 2, A H Ci = {{x,^;y, r])\x2 = = ^2} 

A intersects C2 cleanly in codimension 2, A n C2 = {{x,^;y, r])\xi = = ^1} 

Ci intersects C3 cleanly in codimension 2, Ci n C3 = {{x,^;y,r])\xi = = ^1}. 

C2 intersects C3 cleanly in codimension 2, C2 n C3 = {(x, ^; y, r/)|x2 = = ^2}- 

An(:73 = = C'inC2. □ 

Theorem 5.2. Let F be as in (jl]) with order m. Then F*F can be decomposed as a sum of operators 
belonging in /2"^'°(A, Ci) + /^'"'^(A, C2) + /^"^'^(Ci, C3) + /2"^'°(C2, C3). 

Proof. Recall from Theorem 12.11 that the canonical relation of F drops rank on the union of two sets. Si 
and Ti2- Accordingly, we decompose F into components such that the canonical relation of each component 
is either supported near a subset of the union of these two sets, one of these two sets or away from both 
these sets. More precisely, we let ipi and ip2 be two infinitely differentiable functions defined as follows 
(refer Figure [I]): 

Mx) = I °" Z/'^^''^^,^ • '"f < '\ and Mx) = I °" Z/'^^'''^,^ ■ 
^ ' ^ 1 0, on {{xi,X2) : \x2\ > 2e} ' 1 0, on {(xi,X2) : |xi| > 2e} 

Then we write F = Fq + Fi + F2 + F3 where Fi are given in terms of their kernels 
Kfo = / e-^^a^iV2dw, Kp, = [ e-^'^a^i(l - '(A2)dw, 



Kf^ = J e-^^a(l - Vi)V'2dw, Kp, = J e-^'^a(l - ^i)(l - V2)dw, 

where (p is the phase function of F in (jH). Now we consider F*F, which using the decomposition of F as 
above can be written as 

F*F = F^F + (Fi + F2)*Fo + F^Fi + F^F2 + FIF2 + ^3*^1 + F^F^ + ^2*^3 + F^F (17) 

The theorem now follows from Lemmas 15. 3| 15.41 and Theorem 15.51 below, where we analyze each of the 
compositions above. □ 
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Figure 1. Support of the cutoff functions 'i/'i and -02 • 



Lemma 5.3. Fq, F^F2 and -F^ -^i ^'"^ smoothing operators. 

Proof. We will only prove that Fq and F*F2 are smoothing. The proof for F2F1 is similar to that of F*F2. 
Let ip = ^(/?, where </5 is the phase function in (jj]). 

For 5 = 18e^/h, we analyze Fq according to the following cases: 

(a) {{s, t) : s > 0,0 <t < 2yJ s'^ + - b}. 

For this case, we show that Kp^^ is smoothing. For, on : s > 0,t < 2\/ + /i^ — (5}, ^ is 

bounded away from and hence is a smooth function. Therefore for any m > 

KF,{s,t,x) = Jo:: Mx)Mx)a{s,t,x,u)du;. 

Now by integration by parts, the order of the amplitude can be made smaller than any negative number. 
Therefore Kpg is smoothing. 

(b) {{s,t):s> 0, \t - 2Vs2 + /i2| < 5}. 

For (s,t) in this set, the kernel Kpf^ is identically due to our choice of the function g{s,t) in ([3]). 

(c) {{s,t) : s > 0,t > 2Vs^ + h? + 5}. 

In this case, we have that depending on our choice of the kernel Kp^^ is either identically 
or smoothing. For, if we consider x in the complement of the set (— 2e, 2e)^, then due to the fact 
that supp('(/'i(x)'02(a^)) C [— 2e,2e]^, we have that the kernel is identically 0. Now if we consider 
X G (— 2e,2e)^, then is never vanishing. Then by an integration by parts argument as in Case (|a| 
above, we have that Kp^^ is smoothing. 

Now we consider F2F1. We have 

Kp,{s,t,x)= [ e-*'^^("'*'^)Vi(2;)(l-V'2(2;))a(s,t,x,t^)dw. 
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and 



Due to the cut-off functions -01 and "02 in these kernels, we are only interested in those singularities lying 
above a small neighborhood of the rectangles with vertices (ibe,ibe), (ibe,ib2e), (ib2e,±e), (ib2e,ib2e). 

We have that Kp-^ is smoothing when x values are restricted to a small neighborhood of these rectangles. 
For, as in the previous case, we consider the three cases: For Cases (jaj) and (jcj), the kernel Kp-^ is smoothing 
and the proof is identical as before. For Case ((b]), due to the choice of the function g{s,t), the kernel 
Kp-^ = 0. Therefore F2F1 is smoothing. 

□ 

Lemma 5.4. F^F-^, F2F3 and F^F can be decomposed as a sum of operators belonging to the space 

/2'"(A) + P^^iCi \ A) + /2'"(C2 \ A) + /2'n(C3 \ (Ci U C2)) . 

Proof. Each of these compositions is covered by the transverse intersection calculus. Below we will prove 
for the case of F^F^. For the other operators, the proofs are similar. 

Let us decompose F3 = F3 + Fg + F3 + F3 and = (F/ + F^ + Ff + Ff)*, where the superscripts in 
both these sums denote restriction of F3 and Fi, respectively, to each of the four quadrants. Note that in 
the decomposition of Ff, we stay away from Si by introducing a microfocal cutoff. This is valid because 
the support of the canonical relation of F3 stays away from Si U S2. Then we have 

{Fl)*Fi G I^'^'iA), {FlYFi E l2'"(C7i \ A), {F^yFi G /^'"(Cs \ A) and (F/)*F| g /^'"(Ca \ (Ci U C2)). 

The other compositions can be considered similarly. □ 

We are left with the analysis of the compositions F*Fi and F2F2. This is the content of the next 
theorem: 

Theorem 5.5. Let Fi and F2 be as above. Then 

(a) F^Fi G /2"^'0(A,Ci) + /2™'0(C2,C3)- 

(b) F^F2 G /2"»'0(A,C72) +/2"'°(Ci,C3). 

Proof. We consider F^Fi. The proof for F2F2 is similar. 

We decompose Fi by introducing a smooth cut-off function ^l'3{x) such that ^3(2;) = 1 for xi > e/2 and 
supported on the right-half plane xi > e/4. That is, we write Fi as 

Fi = F+ + Ff, 

where 

F+V{s,t) = j e-''<'^'^^'''''^^i{x){l-i^2{x))i)z{x)a{s,t,x,uj)V{x)dx 

and 

Ff y(s,t) = j e-^^('''*'^'^)0i(x)(l - - iJ:i{x))a{s,t,x,u:)V{x)dx 

Now 

Fi*Fi = (F+)*F+ + (Ff )*F+ + (F+)*Ff + {F^TF^. (18) 

The canonical relation of F^ is a subset of ([8]) with the additional condition that xi > e. Then we have 
that VFF((F+)*F+)' C AuCi. For, we already saw in Proposition O that WF{F*F) C AUC1UC2UC3. 
In our case, imposing the additional restriction that xi > e, the only contributions are in A and Ci. By a 
similar argument, we have that WF{{F^)*F^y C Au Ci. 

Now let us consider the compositions {F^)*F^ and {F^)*F^ . The wavefront sets of these operators are 
of the form (x, ^, y, rj) such that xi and yi have opposite signs. We have already established in Proposition 
15.11 that \xi\ = \yi\ and = |r/j| for i = 1,2. Now with the additional restriction that xi and yi have 
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opposite signs (and therefore and rji have different signs as well), we have contributions contained in 
only C2 and C3. 

The Lagrangian pairs A,Ci and 6*2,6*3 intersect cleanly. Therefore there is a well-defined I^'' class - 
which we will identify shortly - in which each of the summands in (jlSp lie. 

We now show that (Ff )*Ff G P"''^{A,Ci) and that (Ff )*F+, (F+)*Ff G l2™.0(C2, C3). 

We follow the ideas of [5], where the iterated regularity theorem of was used to prove an analogous result. 
The ideas of [5] were recently employed to prove a similar result for a common-offset geometry in [11] . The 
proof we give is similar to the one given in [TT], but the phase function we work with is different. 

We first consider the generator of the ideal of functions that vanish on A U Ci [5] . 

pi = xi-yi, P2 = xl-y2, P3 = Ci-m, Pi = {x2 + y2){i2-m)^ 

Let Pi = QiPi, for 1 < i < 6, where qi,q2 are homogeneous of degree 1 in (^,7/), 93,54 and are 
homogeneous of degree in (^, rj) and qe is homogeneous of degree —1 in (C, ry). Let Pi be pseudodifferential 
operators with principal symbols pi for 1 < i < 6. 

We show in Appendix |A] that each pi can be expressed in the following forms: 

PI = ^^^^^^ds<^ + h2ix,y,s)d^'^ (19) 
P2 = ^^^^^^ds^ + f22ix,y,s)d^^ (20) 

00 

P3 = f3i{x,y,s)ds<^ + ujf32{x,y,s)d^<^ (21) 

Pi = fAi{x,y,s)ds^ + ujfi2{x,y,s)duj^ (22) 

P5 = hi{x,y,s)ds^ + u:h2{x,y,s)d^^ (23) 

P6 = w/6i(x,y,s)a^$-Fa;2/62(2;,y,s)9^$ (24) 

where fij for 1 < i < 6 and j = 1,2 are smooth functions. 

Now the rest of the proof is the same as in [5l Theorem 1.6]. We give it for completeness. 

Let be the kernel of {F^)*F^. This is the kernel in ()10p . but with a there replaced by 
V'i(a;)(l — ilj2[x))ilj'^{x)ilji{y){l — V'2(y))V'3(y)fl(2^) 2/) w). For simplicity, we rename this as a again. 

We then have that a G 5'2™+i ^^^^ 

PiKt{x,y) = I e'''^-'y'''^^a{x,y,s,u)q, (^M^^llAd,^ + f,^{x,y, s)d^^^ dsdu 
ds (e'*(^'?'''*''^)) ^a{x,y,s,u;)fn{x,y,s)dsdu; 
+ Jd^ (^e'*(^.?/.«''^)) ^a{x,y,s,io)h2{x,y,s)dsdio 

By integration by parts 




iH-,v,s,u.)d^ (^a{x,y,s,Lo)fu{x,y,s)) dsdto 
m^,y,s,.)Q^ (^2(x,y,s,a;)/i2(x,y,5)) dsdcol. 



Note that qi is homogeneous of degree 1 in w, and a is a symbol of order 2m + 1, hence each amplitude 
term in the sum above is of order 2m -\- 1. 

Therefore by Definition 13. 7^ we have that PiK^ G Hi°^ for some sq. 



12 AMBARTSOUMIAN, FELEA, KRISHNAN, NOLAN AND QUINTO 

A similar argument works for each of the other five pseudodifferential operators. Hence by Proposition 
13.71 we have that {F^)*F^ G /^''(A, Ci). Because C is a local canonical graph away from S, the transverse 
intersection calculus applies for the composition {F^)*F^ away from S. Hence {F^)*F^ is of order 2m 
on A \ Ci and Ci \ S. Since {F^)*F^ is of order p + / on A \ S and is of order p on Ci \ S, we have that 
p = 2m and / = 0. Therefore {F+)*F+ G l'^"'^'^{A,Ci). Similarly (Ff )*Ff G /2™'0(A,Ci). 

To show that {F^)*F^, {F^)*F{' G /^™'°(C2, C3) we can use the iterated regularity result as above. 

The generators of the ideal of functions that vanish on C2 U C3 are: 

n = xi + 2/1, r2 = ^1 + r]i and p2, pi, Ps, Pe are the same as in (pOj) . (p2|) . (p3|) and (pi]) respectively. 

Four of the functions in the ideal are the same as in the proof above and we can find similar expressions 
for the first two. 

However we will also give an alternate proof below. 

□ 

Proposition 5.6. (Ff )*F+, (F+)*Ff G /^m.o^^;;^, C3). 

Proof. We show for {F^)*F^ . The proof for the other case is similar. Consider the operator R defined as 
follows: 

RV{xi,X2) = V{-Xi,X2). 

This is a Fourier integral operator of order with the canonical relation C2. This is because, 

RV{xi,X2) = j e*("-f) «i?2V^(yi,y2)dydC 
e<-y)<V{-y^,y2)dydi 

ei((xi+j/i)6+{^2-3/2)6)y(y^^y2)dydC 

It is easy to check that canonical relation is C2. 

Now consider the operator F = F^ o R. This is given by 

FV{s, t)= [ e-*'^(^'*'"''^)(^i(l - V2)(l - V3))(-xi, X2))a(s, t, -xi,X2,uj)V{x)dxduj 



Note that (1 — '4>3){—xi,X2) = 'ip3{xi,X2) except in a small neighborhood of the origin. Since ^i(l — ^2) is 
in a neighborhood of the origin, and noting that we can arrange ipi and ip2 to be symmetric with respect 
to xi, we have 

-''^(''*'^''^)(^i(l- V2)V'3))(x)a(s,t,-2;i,X2,w)l^(2;)dxda;. 

Now we have that F*F^ G /^™'''^(A, Ci). In fact the kernel of this operator has the same form as in 
(llOp and the same proof as in Theorem 15.51 applies. Next we use [H Proposition 4.1] to show that 
R*F*F+ G I^'^'°{C2,C3). It is straightforward to check that C2 o A = C2, C2 o Ci = C3 and C2 x A 
(as well as C2 x Ci) intersects T*X x At*x x T*X transversally. Hence the hypotheses of [U Proposition 
4.1] are verified and we conclude that R*F*F+ G /2™'0(C2, C3). Since F* = /?*(Ff )* and {R*f = Id we 
have (Ff G I^^.o^^^, C3). □ 

Since /^'"(A) G /2™'0(A, Ci); /2™(Ci \ A) G /2™'0(A, Ci); /2'"(C72 \ A) G P'^'°{A,C2) and P'^iCa \ 
(Ci UC2)) G /^'"''^(Ci, C3) then Theorem 15.21 follows using Lemmas 15.31 15.41 Theorem 5.5 and Proposition 
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Remark 5.7. Using the properties of the F-' classes, F*F G /2'^'°(A,Ci) imphes that F*F G /2™(A\Ci) 
and F*F G /^'"(Ci \ A). This means that F*F has the same order on both A and Ci which imphes that 
the artifact Ci has the same strength as the initial singularities given by A. Similarly for C2 and C3. Note 
that Ci gives an artifact that is a reflection in the xi axis, C2 gives an artifact that is a reflection in the 
X2 axis, and C3 gives an artifact that is a reflection in the origin. 
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Appendix A. 

Here we give the derivations of ()19p - (j24p . We will work in the coordinate system defined in (|16p . We 
extend the phase function in (jlip to by letting 



$ = ^(yi - s)2 + y2 + (yg -hy + ^ (yi + s)2 + yl + (yg - h) 




Then note that 



(Xi - s)2 +X2 + (X3 - /l)2 + J (xi + s)2 + X2 + (X3 - h) 



dui^\xs=y3=o = du,^ and ds^\x.^=y3=o = ds^. 



A.l. Expression for xi — yi- We obtain an expression for xi — yi in the form 

Ai := xi-yi = M^llAd.q, + /,2(x, y, s)5^$, 

CO 

where fu and /12 are smooth functions. In the coordinate system (|16p . 

Ai = s (cosh p cos (j) — cosh p' cos (j)') 
d^cp = 2s (cosh p' — cosh/)). 



(9,<I> = uj< 



/ cosh p' cos (/)' + 1 cosh p' cos — 1 \ / cosh p cos (p + 1 cosh p cos (p — 1 
V cosh p' + cos (/)' cosh p' — cos / V cosh p + cos (p cosh p — cos 




cosh p' — cosh p' cos^ cp' cosh p — cosh p cos^ ^ 



2uj< 



cosh p' — cos2 (j)' cosh p — cos^ ^ 
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After simplifying, we get, 



2 ( (cosh p — cosh p') (cosh p cosh p' — cos^ (p cos^ (/>') ^ 
I (cosh^ p' — cos^ i;^>')(cosh^ /9 — cos^ (/>) 

(cosh p' cos^ — cosh p cos^ (/)')(cosh p cosh p' — 1) 1 
(cosh^ p' — cos^ (j)'){cosh^ p — cos^ (f)) J 

Now observing that cosh p' — cosh p = and adding and subtracting cosh p cos^ to the second term on 
the right above, we have, 

„~ u (cosh p cosh p' — cos^ (icos^ (i') „ ~ 

H — — duj^ = 

s (cosh^ p' — cos^ (f)'){cosb? p — cos^ 



2uj 



(cosh p' — cosh p) cos^ <j) + cosh p(cos^ — cos^ <?^') | (cosh p cosh p' — 1^ 



(cosh^ p' — cos^ (/)')(cosh^ p — cos^ cf)) 
From this we get 

,, 5s$(cosh^ p' - cos^ (/)')(cosh^p - cos^ (/>) 

cos 0- cos =- — 7 — ; —7 ] -pr (25) 

ILO cosh p(cosh p' cosh p — lj(cos cp + cos q)') 

^di^^ ((cosh p cosh p' — cos^ 0cos^ (j)') — cos^ ^(coshpcosh p' — 1)) 
2ijj cosh p(cosh p' cosh p — 1 ) (cos (j) + cos (^') 

Now note that 

^ s5s$(cosh^ p' — cos^ (/>')(cosh^ p — cos^ (j)) 
^ 2a;(cosh p' cosh p — 1) (cos + cos 

ujdi^<^ ((cosh p cosh p' — cos^ cos^ 0') — cos^ (/)(cosh p cosh p' — 1)) cos (j)' ~ 

H ^ ^ di^^ 

2u} (cosh p' cosh p — 1 ) (cos </> + cos (j)' ) 2 

Now letting = y^, = 0, we see that we have written xi — yi as a combination in terms of and 
as follows: 

s(cosh^ p' — cos^ (/'')(cosh^ p — cos^ (/>) 

2(cosh p' cosh p — 1) (cose/) + cos (/)') uj 

( ((cosh p cosh p' — cos^ </) cos^ (j)') — cos^ (/)(cosh p cosh p' — 1)) cos (j)' ^ ^ 

s (dco'^ 

I 2(cosh p' cosh p — 1) (cos (p + cos (^') 2 J 

We can write the above expression in the Cartesian coordinate system. First, for simplicity, let 



Xi = J{xi - s)2 + xi + /i2 



^2 = Y (2:1 + S)2 + x| + /l2 

with Yi and Y2 being similarly defined with x replaced by y. Then we have 
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2s I \ 2s 



yi 



Yi+Y2 \ ( Xi+X2 \ A( I yi \ Y1 + Y2 

2s M 2s J ^ \ x-^+X2 Y1+Y2 



2 2 



A. 2. Expression for x| — 2/2* Now we write — in form 



A2 := xl - yl = M^lMa^$ + /22(x, (26) 

bJ 



wliere /21 and /22 are smootli functions. A2 in the coordinate system (jl6p is 

^2 = (sinh^ p sin^ (j) cos^ — sinli^ p' s\t? (j)' cos^ 0') 

= (sinh^ p sin^ (j) - sinh^ p' sin^ (p') + (27) 

(sinh^ p' sin^ 0' sin^ 9' - sinh^ p sin^ (j) sin^ 6*) . (28) 

For X3 = 7/3 = 0, (j28p is 0. Therefore we focus only on the term (|27p . which we still denote as yl2) and 
obtain an expression of the form (j26p for this term. 

Using the formulas sinh^ p = cosh^ P ~ 1) sin^ = 1 — cos^ 0, and simplifying, we have 

A2 = ((cosh^ p — cosh^ p) s\T? cj) ~ (cos^ (j) — cos^ (j)') sinh^ p') 

= ((cosh p — cosh p') (cosh p + cosh p') — (cos — cos (j)'){cos cj) + cos 0') sinh^ p') . 

Recall that cosh p — cos p' = and using the expression for cos — cos (j)' in ()25p . and setting 2:3 = ^3 = 0, 
we see that — c^'^ be written in the form ([26j) . 

A. 3. Expression for ^1 — 771. Note that ^1 = dx^^ and ryi = — and so A3 := Ci ~ ^1 = dx^^ + 5^^*^. 
We have 



\y-lT{s)\ |y-7i?(s)K Ma;-7T(s)| |x-7ij,(s) 
In the coordinate system (I16p this is 



2a; 



(sinh^ p' cos (j)' sinh^ p cos \ 

cosh^ p' — cos^ (j)' cosh^ p — cos^ y 



Simplifying this, we get 

(cos (f> — cos 4>'){— sinh^ p cosh^ p' — sinh^ p' cos (/> cos (j)') 



2oj 



+ 



(cosh^ p — cos^ (/>)(cosh^ p' — cos^ (f>') 
(cosh p — cosh p') (cosh p + coshp') cos (f)'{l — cos (/)cos (/>') 



(cosh^ p — cos^ (/>)(cosh^ p' — cos^ 0') 

we can write A^ in the form (|2ip . 



Now noting that coshp — cosh p' = and using the formula ([25]) for cos i;^ — cos (p' and setting X3 = 2/3 = 0, 
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A. 4. Expression for (x2 — 1/2) {S,2 + Using the coordinate system (|16p . we can write A4 := {x2 + 
2/2) (C2 - m) (up to a negative sign) as 



/ "^2 

{x2 - y2)(6 + m) = ^{x2 - 2/2) 1 r 



+ ^2 _^ y2 _^ y2 



\x - ir\ \y - 7t| \y - ir\ 



x\ cosh p 


cosh p' 


cosh^ p — cos^ 6* 


cosh^ p' — cos^ 0' 


X2y2 cosh p' 


2:22/2 cosh p \ 


cosh^ — cos^ 6' 


cosh^ p — cos^ 6 j 


X2 cosh p 


x\ cosh p' 



^ cosh^ p — cos^ 9 cosh^ p' — cos^ 9' 
, 9 9, coshp' 

+ {xi - yi) 9 — ^ 

cosh p' — cos^ ^' 

_l_ 2:22/2 cosh p' 2;2?/2 cosh p 

cosh^ p' — cos^ 9' cosh^ p — cos^ 9 



cosh. 

Here we have added and subtracted ^^^^ pi_^^^i qi ™ the previous equation. Simphfying this we get, 



2uj I 2 

{x2 - 2/2) (6 + m) =— I \.X2 - 2:22/2) 



+ 



(cosh p cosh p' + cos^ 0) (cosh p' — cosh p) 
(cosh^ p — cos^ 0)(cosh^ p' — cos^ 6*') 
cosh p(cos 9 + cos 6*') (cos 9 — cos 0') 
(cosh^ p — cos^ 0)(cosh^ p' — cos^ 0') 



I 2^ COshp^ \ 

cosh p' — cos^ 9' I 

Now note that cosh p' — cosh p = and we already have expressions for cos 9 — cos 9' (Equation ()25p ) 
and for X2 — 2/2 involving combinations of 9^^$ and 

Hence we can write (2:2 — 2/2)(C2 + ??2) in the form of (|2'2|) . Note that our calculation in this section shows 
that 

cosh p cosh p' (cosh p cosh p' + cos^ 9) (cosh p' — cosh p) 



cosh^p — cos^^ cosh^ p' — cos^ 0' (cosh^ p — cos^ 0)(cosh^ p' — cos^ 0') ^^g^ 

^ cosh p(cos + cos 9') (cos — cos 9') 
(cosh^ p — cos^ 6')(cosh^ p' — cos^ 9') 

This will be useful in the derivation of (j24p in lA. 61 below. 



A. 5. Expression for {x2+y2){^2 — 112) • This is very similar to the derivation of the expression we obtained 
for {x2 - 2/2) (6 +m)- 



A. 6. Expression for £,2 — vi- have 
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2 ^2 



^2 - % = ^ 1 r + 1 r - 1 r + i 

cosh^ /9 2 cosh^ p' 



(cosh p — cos^ (cosh p' — cos^ 0')^ 

cosh^ p cosh^ p' 

(cosh^ p — cos^ 9)'^ (cosh^ p' — cos^ 9')'^ 



2 2n cosh^p' 

+ - y2) 



(cosh p' — cos^ ^')^ J 

Now usmg the computations for — 1/2 ^'^d (^2 — 2/2)(^2 + ?/2)5 hi particular ([29]) . we can write ^2 ~ ^2 i'^ 
the form 

^2 - ??i = a;/6i(a;,y,s)9s$ + a;^/62(x,y,s)9^^> 

for smooth functions /gi, /62- 



Appendix B. 

Here we explain the reason for setting g{s,t) = for \t — 2^/s^~+l]?\ < 20e^//i in ([3]). 

In the proof of Theorem 15.21 ~ more precisely Lemma 15.31 - recall that we consider four squares with 
vertices (ibe,ibe), (ibe, ib2e), (ib2e,±e), (ib2e,ib2e). The motivation to choose 5 = as above comes from 
the fact that we want the amplitude term a of F to be for those {s,t) such that the ellipse defined by it 
is contained in a small neighborhood containing these squares. 

One way to find this is as follows: 

Given {s,t), the ellipse (xi — s)^ + X2 + h? + \/ (xi can be written in the form 

(4*2 - lQs')xl + At^xl =t'- At\s^ + h^). 

Note that for this ellipse, the length of the semi-minor axis is always smaller that the length of the semi- 
major axis. The point (2e, 2e) is 2^/2e away from the origin. Therefore let us choose a t for which the 
ellipse passes through the point (0, 3e). The time t is such that 

Hence 

t - 2Vs^ + = 36eV(t + 2Vs^ + h^). 

Since t > and s > 0, we have 

t - 2 + /i2 < 18eV/i- 

This explains the factor 18 in Lemma [5. 31 Now choosing 20 (any number bigger than 18 would do) explains 
our choice of the constant in ([3]). 
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